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We challenge the widely held belief that the cosmological principle is an obvious consequence of 
the observed isotropy of the cosmic microwave background radiation (CMB), combined with the 
Copernican principle. We perforin a detailed analysis of a class of inhomogeneous perfect fluid 
cosmologies admitting an isotropic radiation field, with a view to assessing their viability as models 
of the real universe. These spacetimes are distinguished from FLRW universes by the presence of 
f"^ ■ inhomogeneous pressure, which results in an acceleration of the fluid (fundamental observers). We 

examine their physical, geometrical and observational characteristics for all observer positions in the 
C^l , spacetimes. To this end, we derive exact, analytic expressions for the distance-redshift relations and 

■ anisotropics for any observer, and compare their predictions with available observational constraints. 

As far as the authors are aware, this work represents the first exact analysis of the observational 
properties of an inhomogeneous cosmological model for all observer positions. Considerable attention 
is devoted to the anisotropy in the CMB. The difficulty of defining the surface of last scattering in 
' exact, inhomogenous cosmological models is discussed; several alternative practical definitions are 

presented, and one of these is used to estimate the CMB anisotropy for any model. The isotropy 
constraints derived from 'local' observations (redshift < 1) are also considered, qualitatively. A 
crucial aspect of this work is the application of the Copernican principle: for a specific model to be 
acceptable we demand that it must be consistent with current observational constraints (especially 
anisotropy constraints) for all observer locations. The most important results of the paper are 
^s^j , presented as exclusion plots in the 2-D parameter space of the models. We show that there is a 

region of parameter space not ruled out by the constraints we consider and containing models that 
are significantly inhomogeneous. It follows immediately from this that the cosmological principle 
ON ■ cannot be assumed to hold on the basis of present observational constraints. 

ON 
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I. INTRODUCTION 



In a previous paper Q we found the complete class of irrotational perfect fluid cosmologies admitting an isotropic 
radiation field for every observer, and showed that they form a subclass of the Stephani family of spacetimes that 
includes the FLRW models, but, more importantly, also includes inhomogeneous cosmologies. Here we study the 
' characteristics of these inhomogeneous spacetimes from all observer positions in order to assess their consistency with 
observations, especially with regard to isotropy constraints. In particular, we consider whether the FLRW models are 
the only viable candidates for a cosmological model. A pivotal part of this analysis is the cosmological principle (CP) 
and its relation to the Copernican principle and the observed isotropy of the universe. 

Although it is an assumption based purely on philosophy, the Copernican principle is intuitively very appealing 
[^|,^) : in order to reject an inhomogeneous cosmological model on the basis of its conflicting with the observed isotropy 
of the universe it is necessary to consider all observer positions and to show that for most observers in that spacetime 
the anisotropy observed is too large to be compatible with observations. Actually, in this paper we will, for simplicity, 
require consistency with observations for all observers - our results will thus be rather stronger than is strictly required 
by the Copernican Principle. 

Having adopted the Copernican principle, the question then arises as to whether the observed isotropy of the 
universe, when required to hold at every point, forces homogeneity, thus validating the CP. Well, the nearby universe 
is distinctly lumpy, so it would be difficult to claim there is isotropy on that basis. However, the CMB is isotropic to 
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one part in 10 3 . Together with the Ehlers, Geren and Sachs (EGS) theorem or rather, the almost EGS theorem 
of Stoeger, Maartens, and Ellis [|| this allows us to say that within our past lightcone the universe is almost FLRW 
(i.e., almost homogeneous and isotropic), provided the fundamental observers in the universe follow geodesies (that 
is, as long as the fundamental fluid is dust). What the almost EGS theorem achieves is to provide support for the 
CP without the need to assume the (near) isotropy of every aspect of the spacetime: isotropy of the CMB alone is 
enough to ensure the validity of the CP (for geodesic observers). Since, however, the Stephani models of JlJ do not 
satisfy the conditions of the almost EGS theorem (the acceleration is non-zero), its conclusions do not apply; indeed, 
we know from jjj that there are inhomogeneous spacetimes admitting an isotropic CMB (see also Q). 

A number of inhomogeneous or anisotropic cosmological models have been studied in relation to the CP. The 
homogeneous but anisotropic Bianchi and Kantowki-Sachs models (see also Q §6, and references therein) have 
been investigated with regard to the time evolution of the anisotropy. It can be shown, for example, that there exist 
Bianchi models that spend a significant phase of their evolution in a near-FLRW state, even though at early and 
late times they may be highly anisotropic (again, see §6 of Q and references therein). Most interestingly for the 
work presented here, ^ show that there exist Bianchi models having an almost isotropic CMB but which are actually 
highly anisotropic (and have a large Weyl curvature, emphasising that they are certainly not FLRW), undermining 
the almost EGS theorem of ^ . 

Of the inhomogeneous models that arise in cosmological applications, by far the most common are the Lemaitre- 
Tolman-Bondi dust spacetime s [ L0H12[ ■ These are used both as global inhomogeneous cosmologies - probably the 
most important papers being |l3jQ, studying geometrical aspects, and |l5|-p0t| investigating observational aspects - 
and also as models of local, nonlinear perturbations (over- or under-densities) in an FLRW background See 
also p5{ for a review. Non-central locations are rarely considered in the literature in the analysis of inhomogeneous 
cosmologies, owing to the mathematical difficulties that this usually entails, but |^6| made a study of Tolman-Bondi 
spacetimes from a non-central location and applied their results to a 'Great Attractor' model. Most other non-central 
analyses, though, look only at perturbations of standard FLRW models. 

However, there has been some consideration of Stephani solutions ]27|-]3l|] ) . These are the most general conformally 
fiat, perfect fluid solutions - and obviously therefore contain the FLRW models. They differ from FLRW models in 
general because they have inhomogeneous pressure, which leads to acceleration of the fundamental observers. |52| fitted 
a certain subclass of these models to the first type la supernova (SNIa) data of j}3| using a low-order series expansion 
of the magnitude-redshift relation for central observers derived in |34j, and found that they were significantly older 
than the FLRW models that fit that data. In (3^] (see also ||) we extend those results, using exact distance- redshift 
relations, and show that it would always be possible to find an acceptable fit to any data that could also be fit by 
an FLRW model (for plausible ranges of the FLRW parameters Ho, flo and ^a). The best-fit models were again 
consistently at least 1 — 4 Gyr older than their FLRW counterparts, with models that provide the best fits to the 
newer SNIa data f36f giving an even greater age difference. 

One feature of the Stephani models that has been the subject of much debate is their matter content. The usual 
perfect-fluid interpretation precludes the existence of a barotropic equation of state in general (because the density 
is homogeneous but the pressure is not), although a significant subclass - those with sufficient symmetry - can be 
provided with a strict thermodynamic scheme [p] |3^j3q |. Individual fluid elements can behave in a rather exotic 
manner, having negative pressure, for example (cf. Sec. [ID ). For these reasons, amongst others, [|39f has claimed 



that Stephani models are not a viable description of the universe, but 31, p. 170], argues rather vigorously that this 



conclusion is incorrect, as do we. Cosmological models are often ruled out a priori because the matter fails to obey 
some or all of the energy conditions [ Eopj] ]. However, it has become increasingly difficult to avoid the conclusion 
that the expansion of the universe is accelerating, with the type la supernovae data of ]3(| being the latest and 
strongest evidence for this (see also j42||l3]]). This suggests that there is some kind of 'negative pressure' driving the 
expansion of the universe. For FLRW models, this must correspond to an inflationary scenario, with cosmological 
constant A > 0, or a matter content of the universe which is mostly scalar field ('quintessence' - see [p4[-^7[; see also 
p8| for a discussion of the dynam ical effects associated with A). Such matter inevitably fails to satisfy the strong 
energy condition. We show in Sec. [I D that the subclass of Stephani models we consider contains some which satisfy 
all three energy conditions, and others which do not. In particular, we suggest in Sec. ^ that the models that would 
provide the best fit to the supernova data are models which also break the strong energy condition. 

The main aim of this paper is to discuss the observational effects that arise in spacetimes admitting an isotropic 
radiation field when the observer is at a non-central location (the models we consider have spherical symmetry). 
Principally, we are interested in the anisotropy that this will introduce in observed quantities such as redshift, and 
we wish to compare these anisotropics with presently available observational constraints to see whether it is possible 
to rule out the inhomogeneous Stephani models on the basis of their anisotropy. We also examine the anisotropy 
of the CMB in these models: although they admit an isotropic radiation field, which is usually implicitly identified 
with the CMB, it cannot be assumed that decoupling will actually produce this radiation field in an inhomogeneous 
universe. We therefore suggest a number of ways to estimate the effect of inhomogeneity on decoupling, and determine 
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the resulting CMB anisotropy. Of course, we can always make any anisotropy as small as we like by assuming the 
observer is close to the centre of symmetry. But this is a very special position, so that such a resolution to the 
anisotropy problem would be in conflict with the Copernican principle (although it may be possible to circumvent 
this by invoking anthropic arguments |49| ). In addition to the anisotropics, though, we also consider the constraints 
imposed by measurements of the value of the Hubble constant and t he age of the universe, as well as constraints on 
the grosser features of the distance-redshift relations (see Sec. [II C| ) . While we do not perform fits to the available 
magnitude-redshift data here, we assume the results of |j 35 1, which show that it is possible to fit the observed SNIa 
data adequately with the models we consider. 

In the following section we describe the spherically symmetric Stephani models and discuss their physical and 
geometrical properties, and we present in some detail the particular two-parameter subclass of the Stephani models 
that we will be studying. It will be shown that these models do not have particle horizons, in contrast to standard 
FLRW models. The energy conditions will be used to constrain the model parameters, leaving a manageable parameter 
set to be investigated further. Then, in Sec. |II B| , the transformation to coordinates centred on any observer will be 
derived. In Sec. [II the various distance-redshift relations are presented for the models we consider, and observations 
are applied to constrain the values of the model parameters. The constraints we address are the va lue of the Hubble 
constant, age, the 'size' of the spatial sections (the meaning of which will be explained in Sec. [II C) and, most 
importantly, the anisotropy of the CMB. We show that after applying these constraints there remains a region of 
parameter space containing models consistent with all of the co nstra ints. Furthermore, we demonstrate in Sec. IV 
that many of the models not excluded by the constraints of Sec 



III are distinctly inhomogeneous 
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In Sec. 



Ill E| we 

1) observations. In 



examine the constraints on the acceleration and the inhomogeneity provided by 'local' (i.e., z 
two appendices we derive important results used in Sec. [n| 

Readers with families or short attention spans may wish to skip directly to Sec. VI, where the results are summarised 
and discussed. 



II. THE SPHERICALLY SYMMETRIC STEPHANI MODELS 



Stephani models are the most general conformally flat, expanding, perfect fluid spacetimes. They have vanishing 
shear and rotation, but non-zero acceleration and expansion. In []l| we showed that the irrotational perfect fluid 
spacetimes admitting an isotropic radiation field are a subclass of the Stephani models depending essentially on three 
free parameters and one free function of time. Although the general Stephani model has no symmetry at all, we only 
consider the class possessing spherical symmetry (c = 0, or xo = 0, in the notation of [Q - a full analysis of the 
models of [Q without spherical symmetry is given in Q). The metric in comoving coordinates, from [Q, is 

= (1 + jAr 2 ) 2 J R(t) 2 {dr 2 +r 2 dn 2 ) \ (1) 

V(r,t) 2 \ (l + iAr2)2 ^ ) y K 1 

where c is the speed of light, dSl 2 = d9 2 +sin 2 9d(f) 2 is the usual angular part of the metric and the function V = V(r, t) 
is defined by 

V(r 7 t) = l + ± K (t)r 2 : (2) 

n{t) = l-R. t {t) 2 /c 2 . (3) 

R(t) is the scale factor, and V(r, t) is a generalisation of the FLRW spatial curvature factor in isotropic coordinates 
(note that in |j| R(t) is included in V). Since k is a function of t the spatial curvature can vary from one spatial 
section to the next. In fact, it is possible for a closed universe to evolve into an open universe, or vice versa, in stark 
contrast to FLRW models 

We have already restricted ourselves to spherically Stephani models, but we need to reduce the parameter space 
further by introducing a form for R(t) that depends on just a few parameters. To this end, we limit attention to the 
two-parameter family derived in Sec. IV of |5lf| . Thus we choose 

R(t) = ct(at + b) (4) 

where a and b are the free parameters and 

A^l-b 2 , (5) 
so that (as can easily be seen by direct calculation) 
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K(t) = A R(t). 



(6) 



We will henceforth refer to these models as Dabrowski models. We require that after the big bang (and before any 
big crunch) R > 0, which forces 

b > (7) 

In (f|)-(||) we retain factors of the speed of light, c, to facilitate comparison with the references given above 
and with observations. The units we will use are as follows: [c] = km s -1 , r is dimcnsionlcss, R is in Mpc and 
[t] = Mpc s km -1 = [1/Hq], so that Ja] =km s -1 Mpc = [Hq] and b is dimensionless. Note that these units 
are slightly different to those used in [[34| because the parameters a and b in that paper contain a factor of c (so 
[b] = [c], etc.). This explains the appearance of c in (|J). 

We will use T to denote the coordinate time of a specific epoch of observation along some observer's worldline 
(i.e., the coordinate age of the universe), again in Mpc s km" 1 , and r to denote proper time along a particular flow 
line. When we state ages they will generally be given in Gyr: TjGyr] as 978T[Mpc s km -1 ]. 

Formulae for the expansion and acceleration may be found in M , and will be introduced as required, as will formulae 
for the energy density and pressure. 



A. Geometry 

The metric (|l]) is manifestly conformal to a Robertson- Walker metric with curvature A. If we multiply through 
by the conformal factor we see that the spatial sections really are homogeneous and isotropic, but the actual spatial 
curvature is time-dependent and is given by the curvature factor «(t) in V(r,t). Thus, at any time t whether the 
spatial sections are closed, open or flat depends on the sign of k{£). If, at some point during the evolution of the 
universe, t = — (b ± l)/2a, then the curvature changes sign, as can easily be seen from (||) and (^). This does not 
happen in FLRW models, where the spatial curvature, k, is fixed. The distinction between the time-dependent true 
geometry (k) and the fixed conformal geometry (A) should be borne in mind throughout what follows. 

The metric (|l|) is not in its most advantageous form. The conformal geometry of the models is most easily studied 
by changing from the stereographic coordinate, r, to the 'angle' coordinate, x ( see equation (5.15) of |4lj), appropriate 
to the value of A. Furthermore, for models with closed spatial sections (which will be our principal concern here) it 
is more convenient to choose a radial coordinate that is better able to reflect the fact that light rays can circle the 
universe many times. In such models the spatial surfaces have two centres of symmetry, r = and r = oo. Physically, 
there is nothing extraordinary about the point r = oo: it is not infinitely far away from the centre, and it is quite 
possible for light rays to pass through it. This last point is particularly important for subsequent discussions, so we 
make the coordinate change (valid when A > 0, although see below) 



2 
A 



X 

tan — 
2 



(8) 



Then r — > oo as \ ~ * 71 • As a coordinate \ is restricted to the range < x < However, it will prove convenient to 
use x not j us f as a coordinate but as a parameter along light rays. In the latter role its value can increase without 
bound as the rays circle the universe many times. Strictly speaking, we should distinguish these two uses, but it 
should not lead to confusion. The absolute value is taken in (p|) so that, when x increases beyond n, r remains 
positive. Using (||) in the metric (Q) gives 



^ = WfcW {~ c2dt2 + ^( d X 2 +^xd^) \ . CM 



where (using (|| and |^) 



W( X ,t) = cos 2 f V(r(x),f) = cos 2 | + ^sin 2 f = 1 - ^ sin 2 f . (10) 

The conformal factor (1/W) is non-singular for all x if a < 0. When a > singularities W = correspond to spatial 
and temporal infinity, and indicate that the universe has 'opened up'. As the universe opens up and the sections 
become hyperbolic the coordinate x represents a conformal mapping from a hyperbolic surface onto a sphere: spatial 
infinity will then correspond to some finite value of x < tt- For a more detailed explanation of this see Theorems 4.1-4.4 

of a 
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We can easily calculate the acceleration in these coordinates. It has only a radial (x) component: 



2 W,x 2ac R . 

% = -C 2 -^ = — -SHI*. (11) 

A simple calculation shows that the acceleration scalar, which we will need below, is just 

,i^(u a ^) 1/2 = ^sin X (12) 

and is therefore time-independent (which, it turns out, is a necessary condition for a shear-free, irrotational perfect 
fluid to admit an isotropic radiation field - see Note that the units of u are: [it] = [c/t] = [cHq]. 

For completeness we mention that when the Dabrowski model is conformal to an FLRW spacetime with hyperbolic 
geometry (A < 0) the coordinate transformation is obtained from that just given by replacing trigonometric functions 
with their hyperbolic equivalent. 



B. Non-Central Observers 



Since the purpose of this paper is to provide insight into the observational characteristics of the Dabrowski models 
for all observer positions, we must find expressions for the distance-redshift relations and other observable properties 
of the models from any point. For a general inhomogeneous metric this is far from trivial, but the Dabrowski models 
have features that make this tractable (rather simple, actually). In particular, they are conformally flat. As has 
already been noted, in equation (^|) the part of the metric in braces is exactly the form of the FLRW metric in 'angle' 
coordinates, so that the Dabrowski models are manifestly conformal to the (homogeneous) FLRW spacetimes. This 
means that there is a group of transformations acting transitively on surfaces of constant time that preserve the form 
of the FLRW part of the metric (but not the conformal factor). If the spatial sections are closed, flat or open (in 
the conformal sense, i.e., according to the value of A), the transformations are rotations, translations or 'Lorentz 
transformations', respectively. After such a transformation the metric will have the form of an FLRW metric centred 
on the new point, multiplied by a modified conformal factor. 

As will be shown in Sec. [II D| , we will be dealing exclusively with closed models (A > 0), and so will concentrate 
on this case. To find the coordinate transformation to a non-central position, we perform a rotation of the spatial 
part of the metric, moving the origin (x = 0) to the point x = "0 (V 7 is the observer's position in what fo llow s) . In 
appendix [A| we derive this transformation. The old x is given in terms of the new (primed) coordinates by ( |A6| ) . The 
conformal factor ([To]) then becomes (dropping the primes on x and 9) 

W — > W(x, ip, 0;t) = 1 — (1 — cos -0 cos x + sin ip sin x cos 9), (13) 

cA 

while the rest of the metric retains its original form (but now in terms of the new coordinates). This transformation 
makes the study of our inhomogeneous models significantly easier, and allows us to find exact observational relations 
valid for any observer. 



C. Lookback Time and the Horizon 



If we calculate the lookback time in our models (i.e., the time, t, at which a galaxy at some position y emits the 
light that the observer sees now at time T), which we can do directly from the metric ([)]) - see appendix [b|- we get, 
for any observer, 

t(x) = — 7= ■ (14) 

(aT + h) exp (bx/y/A) - aT 

(Note that this function is continuous through x = tt-) Now, t — > if and only if x — ¥ oo. Thus, the whole of 
the big-bang surface is contained within the causal past of every observer in the spacetime and there is no horizon 
problem for the Dabrowski models. This is in sharp contrast to FLRW universes: at early times the particle horizon 
is finite and contains only a small part of the big-bang surface, so that widely separated points can share no common 
influences. See Figs. 17 and 21 in plfl. 
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D. Matter Content and Energy Conditions 



The Stephani models do not have an equation of state in the strict sense, with the relationship between pressure 
and energy density being position dependent. Along each flow line, however, there is a relation of the form p — - 
see fl30[ . The particular models we are using have an equation of state at the centre of symmetry (or everywhere 
in the homogeneous limit a — > 0) of the (exotic) form p = — It was shown in Q), though, that these models do 
admit a (more general) thermodynamic interpretation (that is, temperature and entropy can be assigned in such a 
way that they are functions only of pressure and energy density) . The matter content of these models with regard to 
the natural (comoving) velocity field is a perfect fluid with energy density 



8irG 3 



and pressure given by 



R{tf 



■ > H^tfi), (18) 



3^ V V *Gl* A 2 

where ( |l5| ) has been used to express the i-dependence of pressure in terms of the density. In other words, we 
have a position-dependent 'equation of state' of the form p = — \^c? + e(%)/i 1 / 2 . The appearance of — |/x as the 
dominant contribution (at early times, at least, when R is small so that, from (|l5|), /i is large) immediately suggests 
a quintessential or scalar field interpretation of the matter Q-|^J , although it has been shown |piq , |53[ that cosmic 
strings also give rise to this EOS. Although the energy-momentum tensor of the Stephani models has the perfect fluid 
form, the interpretation of the matter as an actual fluid is by no means required; other interpretations may also be 
valid, and the fact that there is no true EOS might even suggest a two-component interpretation. For the moment, 
though, this is as far as we will go to provide a physical motivation for the matter in the Dabrowski models: in this 
paper we are only interested in the observational consequences of the geometry - although when we impose the energy 
conditions below we show that the matter is certainly not obviously unphysical. 

We can see that there are singularities of density and pressure as R(t) — ► (i.e., at t = 0, —b/a), which correspond 
to the big bang and crunch for these models (the metric becomes singular at these points) - see We can also 
have a finite- density singularity, where only the pressure becomes singular. This happens when r — > 2/y/—A. Such 
infinite pressure is clearly not physical, so we can reject models with A < 0. Models with A = are studied in f$2l, 
but it is difficult to compare them directly with A > models due to the different geometries of the spatial sections 
so we will not consider them in this paper, and we are left with A > 0, i.e., b < 1 - see (||). As we explained in Sec. |0] 
the natural choice of positive R(t) after the big bang ensures that b > 0. 

Having calculated the pressure and density of the fluid we can now investigate its physical viability through the 
energy conditions. It is more convenient to use the original stereographic coordinates for this (i.e., the expression ( |l6| ) 
for the pressure), since we wish to consider all values of A, until we find reasons to the contrary. The weak energy 
condition states that [i > and [ic 2 +p > 0, whereas the strong energy condition requires in addition /ic 2 + 3p > 
(see, for example, p0[ for a discussion). The weak energy condition rules out models with finite-density singularities 
(i.e., with A < 0), because such models would contain regions (for r > 2/V— A) where /ic 2 + p < 0, as can easily 



be seen from (16). The weak energy condition also implies that V > 0, but this is always true since V — > only at 
spatial (and temporal) infinity (even though the coordinates themselves may be finite). The strong energy condition, 
however, implies that n(t) > A for all t. From (^) we can see that this is equivalent to a < (since R > 0), so the 
models must have a big crunch (R(t) is an 'upside-down' quadratic). 

The dominant energy condition is more interesting. It states that \p\ < fic 2 , from which ( |l5| ) and ( fl8|) immediately 
give 

1 V(r,t) 

< Y < 1- 

"31 + iAr 2 " 

The left inequality requires only that A > 0, which rules out infinities in the pressure. The inequality on the right 
says that for all t and r 
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(n(t) - 3A)r 2 < 8 (19) 

must hold. This condition is always true for a > (see (g)), as long as A > 0. When a < 0, r is unbounded for A > 
(because K is then positive - see so the left hand side of ( |l9| ) must always be negative, n(t) < 3A. It is easy to 
see that R(t)/c < ~b 2 /4a, so, from ||, n{t) < 1. Then k < 3A = 3(1 - b 2 ) for all t provided 

6 < y| ~ 0.82. (20) 

For b larger than this the dominant energy condition will be broken at some time, in regions of the universe at large r. 
We will not consider further the intricacies of this. A glance at the exclusion diagrams, Figs. p| and |lO|, shows that ( p5| ) 
does not eliminate a significant area of the allowed region. In light of this we will, for the moment, overlook ( p0| ) and 
investigate the properties of all models with < b < 1. 

To summarise: we have used basic physical requirements, such as the occurrence of a big bang and the avoidance of 
pressure singularities, along with the energy conditions to restrict the ranges that the two parameters a and b (or A) 
can take. The results are: 

a<0, < < 1 (i.e., < A < 1) (21) 



(we reject 6=1 for simplicity, as explained above, and we refrain from invoking (20) until Sec. VI) 



III. CONSTRAINING THE MODEL PARAMETERS USING OBSERVATIONS 

So far we have considered only the 'global' physical properties of Dabrowski models, but to really assess their 
potential viability as cosmological models it is necessary to confront them with observations. In this section we 
derive the distance-redshift relations that form the basis of the classical cosmological tests and compare them with 
available observational constraints to see whether any regions of parameter space are capable of providing a fit. We 
will impose constraints on the value of H 0l age, size (the meaning of which will be explained below) and the anisotropy 
of the microwave background, leaving the wealth of data available from galaxy surveys and high-redshift supernovae 
for consideration in a future paper; the complexities involved in interpreting such data and applying it to idealised 
cosmological models require separate treatment. 

Deriving the observational relations (redshift, angular size or area distance, luminosity distance and number counts) 
means relating the coordinates and metric functions to observable quantities. This requires knowledge of the observer's 
motion (4- velocity) , which can, strictly speaking, be specified independently of the background geometry. However, 
the Dabrowski models contain perfect fluid, so we will identify the observer's motion with the fluid velocity. We are 
not obliged to do this, and, given the strange form of the matter, it might be thought advantageous to instead assume 
that observers (i.e., galaxies) constitute a dust-like test fluid moving freely through the spacetime whose geometry is 



determined by the exotic matter. It will become clear in Sec. Ill D that if we were to make this assumption a large 
dipole anisotropy in the CMB would result (although the dipole in Hq would be eliminated - see ( |30| ) with ii = 0) 
because such a flow will, in general, have a significant velocity relative to the Dabrowski fluid flow, which, it will turn 
out, is very nearly in the rest frame of the CMB everywhere. 

First we consider redshift. In general, it is no simple task to find analytic expressions for the redshift in any 
cosmological model; derivations usually rely on symmetries of the spacetime or other simplifying factors to solve the 
equations of null geodesies. Here we take advantage of the conformal flatness of Stephani models (or rather, of the fact 
that the Dabrowski models are manifestly conformal to FLRW spacetimes), although we can also derive the redshift 



formula as a time-dilation effect; these procedures are outlined in appendix H. Using (B6) and (Bj) we find 



where Rq = R(T), Wo = W(ip,T) and t and x are related by equation (|14J). Using ( |13[) this is 

R ( 1 2a . 2a 1 

1 + z = — ~ cA _ C0S ^ C0S X) ~ ^ sin V> sin x cos 0j ; (23) 

showing that, for objects at any fixed x, the inhomogeneity of universe manifests itself in the redshift as a pure dipole 



in angle around the sky (cos 9 term). This will be important in Sec. HID 
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For metrics with spherical symmetry about the observer the angular size (and area) distance is given directly from 
the coefficient in front of the angular part of the metric, because symmetry ensures that for radial rays 9 and 4> are 
constant along the trajectory. For our models we do not have spherical symmetry about every observer, but the metric 
is everywhere conformal to a spherically symmetric metric, as can be seen from (^). Since null rays are not affected 
by the conformal factor they also remain at fixed 6 and 0, so we can again obtain the angular size distance ta from 
the coefficient of the angular part of the metric: 

r ^'*< 0) = rx^n fl J sinxl (24) 
v Aw(ip, Xi 0; t) 



(again, x an d t are related by (|14|) ; the modulus signs around sin x ensure that ( |24|) is valid even when x is treated as 
a parameter along light rays and takes on values > ir - see Fi g. H ). We can, for the first time, find the exact angular 
size-distance relation parametrically by combining equations ( p4|) and (|2^) , which bypasses the power-series method 
of ]54|| , This is valid for all sources seen from any observer position. 
Luminosity distance is related to ta by the reciprocity theorem: 

r L = (l + z) 2 r Al (25) 

see |^,^5|. This then allows the magnitude-redshift relation to be determined in the usual way: the apparent magni- 
tude m of an object of absolute magnitude M is given in terms of the luminosity distance by 

m-M- 25 = 5 log 10 r£. (26) 

The task now is to limit a, b, and T using present observational constraints. A full discussion of each constraint is 
made in the following sections, and it is followed by exclusion diagrams showing the regions of parameter space for 
which a and b give a plausible cosmological model for all observer locations in these models. 



A. The Hubble Constant 



The expansion rate of the universe has been measured with reasonable accuracy. The Hubble constant is believed 
to lie in the range 50 < Hq < 80 km s _1 Mpc -1 , and we will use these limits to constrain the Dabrowski models. The 
Hubble parameter for these models, which is related to the volume expansion 0, is independent of position 0J31I]: 



H _Q_Rt{t) H - RtiT) (97) 

We can use this to place constraints on the time at which observations are made: for our models H decreases 
monotonically, so it will only lie in the observed range of Hq for some range of T. For any observer with coordinate 
age T, we require 

5o<J^<so. m 

Usually, for simplicity we will choose a specific value for Hq (almost invariably that which produces the 'worst case'). 
Then we can solve (|27]) for T. 

However, when Hq is actually measured, it is not necessarily equal to the expansion rate. What is measured in 
practice is the lowest order term in the magnitude-redshift relation, which gives the measured Hubble constant, H™: 



Tjm 2 ^aUb 
Hq =c 



{u c k c f 



(29) 



where k a denotes the wave- vector of the incoming photons |55| . Equivalently, we can consider the gradient of the 
redshift-area distance curve at the observer (cf. p6[). If we measure the magnitude-redshift relation for objects in some 
direction, then, expanding V a i*b in (^9| ) in terms of the kinematical variables rotation io a b-, shear a a b, acceleration u a 
and expansion O and using (|27|), we find (since a a b = for Stephani models) 

n 

H% l (6) = H Q - -cos 0, (30) 
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where it, is the acceleration scalar and 8 is the angle between the acceleration vector and the direction of observation 
(which is opposite to the direction in which the photons are travelling). Since the acceleration is non-zero in the 
Stephani models there will be a dipole moment in H™. The size of this in the Dabrowski models is given directly 
from (|l2] ) . If this is large in any model we can probably reject that model because a large dipole moment in Hq is not 
observed. However, very nearby it is difficult to measure Hq accurately due to peculiar motions and the discreteness 
of galaxies. There is a dipole moment in observations of somewhat more distant objects, which is assumed to be due 
to the fact that the Local Group is falling into the potential well produced by Virgo and the Great Attractor. The 
question is: what upper bound can be placed on the acceleration by observations? This issue was first raised in [[l], 
and will be discussed in more detail in Sec. Ill E . For now we simply use (pTj) to constrain the epoch of observation, T. 



B. The Age of the Universe 

The original inspiration for fl32| to study Stephani models was the potential resolution of the age problem that they 
provided, which at that time seemed to be virtually insurmountable within the framework of FLRW models (even 
when a non-zero cosmological constant was invoked) : the high measured value of Hq suggested an age, tq, of at most 
about 11 Gyr for an FLRW cosmology, whereas globular cluster ages were thought to be up to 12-13 Gyr p2| , ^5| ,p| 
showed that, for the particular Stephani models they considered, this apparent paradox disappears: the Dabrowski 
models have ages that are consistently 1-4 Gyr older than their FLRW counterparts (for an observer at the centre 
of symmetry, at least). However, the age problem has recently been alleviated by a recalibration of the RR Lyrae 
distance scale and globular cluster ages in the light of Hipparcos |39), which has reduced the globular cluster ages 
considerably, to ~ 10 Gyr. The fit is still marginal, but the new ages are generally accepted as they allow a flat FLRW 
model to fit the observations provided that Hq < 67 km s" 1 Mpc^ 1 . 

We certainly require, then, that at the epoch of observation our models are older than 10 Gyr. However, we will 
also consider the stronger constraint tq > 12 Gyr, partly to be conservative, but also because the diagrams for the 
12 Gyr constraint are often clearer. 

The age of the universe according to an observer at position tf> and at coordinate time T is simply the proper time 
elapsed from the big bang (t = 0): 



_ f T dt A f A - 2aTS (bS + VT= b 2 C 2 ) 1 

T ° 7 W(il>, t) ~ A\a\SVl - b 2 C 2 n \ A - 2aTS (bS -VI- b 2 C 2 ) J ( ' 

(for a ^ and if> ^ 0, otherwise W = 1 and r = T), where S = sin(?/>/2), C — cos(ip/2), and we take the value of T 
given by the solution of equation ( p7| ) as our constraint on the coordinate time for any specific Hq. 

In Fig. [IJthe tq = 12 Gyr contours of the proper-age function are plotted for Hq = 60 km s _1 Mpc -1 and for several 
observer positions, ip, showing how the age of an observer varies with ip for different parameters a and b. The shaded 
regions contain models that are less than 12 Gyr old for at least one of the observer positions. This demonstrates that 
the proper-age of an observer is smallest at the antipodal centre of symmetry, ip — tt. Consequently, we will always 
use proper- age at ip = tt to constrain the model parameters. We could weaken this constraint by requiring only that 
most observers are old enough, which would allow us to consider instead the age of observers at ip < n/2 while still 
satisfying the Copernican principle (half of the observers would lie in this region) . For simplicity, though, we will not 
do this here. 

Finally, in Fig. || we show the age exclusion plot for the models (based on proper-age at ip = n). We use three 
values of Hq: 50, 60, and 70 km s _1 Mpc -1 , and a proper-age of 10 Gyr (although the limits for 12 Gyr are also 
indicated). The shaded regions are excluded. It can be seen that unless we require the universe to be particularly old 
or the expansion rate high there is still a significant region of parameter space that cannot be excluded on the basis 
of age. 

It should be noted here that these plots are meaningless for b — 1, because then A = and the model is conformal 
to an FLRW model with flat spatial sections, for which \ 1S n °t a good coordinate - as can easily be seen from (0). 



C. Size and the Distance-Redshift Relation 



When the spatial sections of a cosmological model are closed and there is no horizon problem light rays may circle 
the entire universe, perhaps many times. This is the case for the models we are considering here. What will the 
signature of this be in the various distance-redshift relations? The paths of light rays are determined by the conformal 
geometry of a spacetime, and it can be seen from (g) that our models are conformal to closed, static FLRW spacetimes. 
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FIG. 1. ro = 12 Gyr exclusion plot for different observer positions, for Ho = 60 km s Mpc . In the shaded areas the 
models are not old enough, ip = tt is clearly the most restrictive case. 

It follows that light rays from a point directly opposite the observer (i.e., from the antipode, X — n ) wm spread out 
around the universe isotropically from the antipode until they pass the 'equator' (\ — t/2), where they will begin 
to converge and be focused onto the observer. As a result, a point source positioned exactly at the antipode will fill 
the entire sky when seen by the observer, so that its angular size distance, ta — (physical length/apparent diameter), 
is zero. Similarly, the refocusing of light onto a point produces an infinite flux at the observer, and therefore the 
luminosity distance, rj,, is also zero (m ~ log 10 — -co). It is obvious that whenever the light rays travel through 
a parameter distance x that is an exact multiple of tt, rA = r^ = 0: this is reflected by the factor of sin^ in (pi]). 

This effect can be seen clearly in Figs, ^-j^, where we show the two principal measures of distance as they vary with 
coordinate distance x or redshift. Viewed as a function of x> i n Fig. ||, the zeros of the angular diameter distance 
occur at multiples of tt for all model parameters. Looked at in terms of redshift, though (Fig. [!]), it is clear that 
for small b the zeros are much closer together than for larger b, with the first zero occurring at z « 1 for b = 0.25. 
Figure |5| shows the luminosity distance-redshift relation, for comparison. These effects are not as unusual as they 
look, and can be found also in FLRW geometries for models with positive A - see §4.6.1 in and references therein. 

Can we rule out such apparently aberrant behaviour? Theories of structure formation are fairly well developed 
(see p3| for a thorough discussion and references), and the evolution of galaxies and the star formation rate (SFR), 
while not accurately known, arc at least qualitatively understood. In particular, the SFR, which is very important 
for determining the luminosity of distant, young galaxies, is believed to fall off beyond z ~ 2 | pl"] , p2fl . As a result, 
one could argue that there will be relatively few bright objects beyond some redshift zsf that corresponds to the 
epoch at which galaxies 'turned on' and the SFR began to increase significantly. This would mean that the zeros in 
the distance-redshift relations would be essentially unobservable if they occurred at redshifts larger than zsf because 
there would be no luminous objects to be seen magnified in the sky, whereas if the zeros occurred at lower redshifts 
than zsf one could reasonably argue that there ought to be some signature of this in the observations. Since galaxies 
have only been observed (in the Hubble Deep Field, for example) with redshifts up to z ~ 5, and quasars have only 
been seen out to a similar redshift, we take zsf = 5. 

The constraints imposed by larger zsf can be inferred from Figs. ^| a nd fij . (For example, if zsf ~ 30, the redshift 
at which it is suggested that the very first luminous objects may appear |63| , Fig. ^, gives b > 0.7 for a = —1.) There 
are, of course, a number of factors that rather cloud these arguments, in particular the possible effects of extinction, 
which arc not well known, as well as the fact that the number of objects 'near' to the antipode will be small, because 
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FIG. 2. The age exclusion diagram for various Ho and proper age r = 10 Gyr. The shaded region represents the prohibited 
area. Also shown as dashed lines are the age limits for r = 12 Gyr. (Note that the region excluded for Ho = 70 km s _1 Mpc -1 
contains the excluded regions for lower Ho - the progressively darker shading indicates this.) 



b=0.7 




FIG. 3. Area distance from the centre as a function of \ for two values of b. T — 15 Gyr, a = — 1. 
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FIG. 4. Area distance from the centre as a function of redshift for the same parameters as Fig. |§| For small b the angular 
size distance oscillates far too rapidly. 



b=0.7 




FIG. 5. Luminosity distance from the centre as a function of redshift for the same parameters as Figs. ^ and | 



12 



b 




the spatial volume is diminished as a result of the geometry. Interestingly, number counts in the HDF show a dip 



at z ~ 2 [64 6Jj] which would occur in any model for which the antipode is at this redshift. However, this dip is 
generally believed to be spurious fl66| . 

If these arguments are not completely convincing, then at a simpler level the fact that the observed magnitude- 
redshift relation is known accurately out to z ~ 1 from type la supernovae [Sfl, and is certainly not dipping down, 
allows us to say that there is no zero of luminosity distance below z ~ 2, say. We therefore also consider the constraint 
that results from requiring that there are no zeros below 2 = 2. 

We wish, then, to constrain the parameters of our models by rejecting any models for which the first zero in the 
distance- redshift relations occurs at z < z w , where z^ = 2 or z w = Zsf = 5. Using (^2|), this means 

1 + Z{X = 7T) = W0 R(U) >1 + x " (32) 

where R — R(T), Wo = W(tp,T) (the conformal factor at the observer) and t„ denotes the lookback time ( |l4| ) 
at x = 7r. Again we determine the epoch of observation (i.e., the observer's coordinate time T) using (^7|). The 
solution of ( |32] ) for a and b is shown in Fig. ^ as an exclusion diagram. The effect of this constraint is to rule out 
small values of b, for any a. This is a reflection of the fact that, loosely speaking, b measures the 'size' of the universe: 
at early times the scale factor goes as bt, so that when b is small the spatial sections are small, light rays don't take 
long to travel from antipode to observer, the scale factor changes relatively little during this time and the redshift of 
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the antipode (which is dominated by Ro/R^) as in FLRW models) is small. 

As a coda to this section we consider the effect of demanding that the first zero of ta is effectively unobservable as 
a result of being 'hidden' behind the CMB. Figure ^| shows how the redshift of the first zero of Ta{z) varies with b. If, 
instead of choosing = zsf as our primary constraint, we want the first zero of Ta{z) to happen at a redshift large 
enough for the universe to be opaque (i.e., before decoupling), then Fig. ^ shows that b must be quite close to unity. 
This figure also allows the extent to which values of b are excluded for any z v to be estimated. 



D. The Microwave Background Anisotropy 

The CMB is observed today to be a blackbody at a temperature of To = 2.734 ± 0.01K, with a dipole moment of 
Ti = 3.343 ± 0.016 x 10~ 3 K and quadrupole moment as large as T2 — 2.8 x 10 _5 K (see |67]] for details and references). 
It was emitted at a time when the radiation was no longer hot enough to keep Hydrogen ionised, causing it to decouple 
from matter, which happens at T^ec ~ 3000K. Idealised cosmological models do not have realistic thermodynamics 
(that is, they do not, in general, accurately describe the thermodynamic evolution of the gas and radiation mixture 
that fills the real universe). In FLRW models the epoch at which decoupling occurs is simply defined to be that 
corresponding to the redshift necessary to shift the temperature at decoupling to the observed mean temperature 
of the CMB, To. From the redshift relation applied to the temperature of a blackbody (T will be used to denote 
temperature in this section), 

T obs = ^L, (33) 
1 + z 

we infer that the CMB is formed at a redshift z w 1000. This definition is fine for homogeneous models, leading to a 
consistent definition of the time of decoupling for every observer at the same cosmic time, but raises an interesting 
point for the inhomogeneous Dabrowski models, because the redshift depends on both the observer's position, ip, and 
the angle around the sky, 9. If we simply define the redshift of the CMB at any point to satisfy ( ]33| ) with T Q b s = T 
then, by definition, we obtain a perfectly isotropic CMB for that observer, but we must choose a different emitting 
surface for each different observer. Such an observer-based definition of the CMB surface is clearly unsatisfactory. 
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In fact, we know already from the results of |jj that the models we consider admit an isotropic radiation field for every 
observer, which is implicitly identified with the CMB. However, this identification overlooks any physics underlying 
the production of the CMB. When, as with the original EGS theorem, the models under consideration turn out to 
be FLRW this is acceptable, since it can be assumed that the homogeneity applies to the production of the CMB: at 
some moment of cosmic time decoupling occurs everywhere throughout the universe. Unfortunately, when the models 
admitting an isotropic radiation field are inhomogcncous this is no longer acceptable, and some consideration must 
be given to the production of the CMB and how this may be affected by the different conditions at different places in 
the universe. To give detailed consideration to the physics of decoupling in the non-standard Dabrowski cosmologies 
would take us beyond the scope of this paper, so instead we consider several alternative, pragmatic definitions of what 
is meant by the CMB surface, which, while not fully capturing the physics of decoupling, at least allow the influence 
of inhomogeneity at the time of decoupling to be estimated (it is important in these definitions to disti nguis h between 



the dominant exotic matter that is responsible for the geometry of the Dabrowski models - see Sec. II D - and the 
putative 'real' baryonic gas that decouples): 

1. If we avoid all consideration of the physics of decoupling, we could simply assume that it happens at such 
an early time that we could define the CMB to be free-streaming radiation 'emitted at the big bang', as is 
effectively assumed in the EGS and almost EGS theorems QJI|. This only really makes sense if there is some 
natural definition of the radiation field at early times. For example, if the model is homogeneous and isotropic 
at early times, we can define a homogeneous and isotropic radiation field. Our models have exactly this property 
of homogeneity at early times (as can be seen from ([lO]), W — > 1 as t — > 0). For our models this definition will 
result in a perfectly isotropic CMB for every observer jlj] (although its observed temperature will be position 
dependent); 

2. Ideally we would like to define decoupling in terms of the thermodynamics of the baryonic gas (at a fixed 
temperature T ~ p/p, say). However, the Dabrowski matter is not an ideal gas and its pressure and density are 
not that of the real baryonic gas, so the utility of this definition here is limited; 

3. In general inhomogeneous models we could choose a fixed value of some physical quantity such as density which 
would allow us to estimate the degree of inhomogeneity - for the Dabrowski models this is equivalent to || 
because the density is homogeneous on cosmic time surfaces; 

4. We could choose a surface of constant cosmic time, t — t CMB - since the Dabrowski models possess a cosmic 
time coordinate with respect to which only the pressure is inhomogeneous this is a natural extension of the 
FLRW definition; 

5. We could take a surface of constant proper time (based on the assumption of some common evolution for the 
ideal gas component at different positions) - see equation (|3l|). 

We will not consider |l| here for the reasons outlined above. Also, the homogeneity of Dabrowski models at early times 
means that for small t the proper-age is virtually identical to the coordinate time (equation ( |3l| ) with W ~ 1). It turns 
out that for times of observation that reproduce the observed Hq any reasonable definition of the CMB surface puts it 
at an early time, which means definition || is virtually identical to We will therefore define the CMB according to ^ 
in this section. This amounts to assuming that the early homogeneity allows for homogeneous physics on the CMB 
surface just as in FLRW models, and that the small inhomogeneity that is present affects only the redshift of points 
on the CMB surface. Note that the anisotropy of the CMB that arises from this definition of the CMB surface does 
not conflict with the results of Q| : a radiation field that is isotropic for every observer may still be defined, but any 
realistic process that gives rise to the CMB must reflect the inhomogeneity of the universe at the time of decoupling; 
the CMB anisotropy we derive in this section results essentially from the (small) inhomogeneity on the CMB surface. 

It still remains, though, to decide exactly which surface of constant cosmic time the CMB originates from. Consider, 
for an observer at position ip, the temperature distribution on the sky that the CMB would have if it were emitted 
from the surface t = t CMB (related by the lookback time formula ( |l4| ) to some distance Xcmb)- 

ToMXcm Bi 0)= Tdcc -. (34) 

1 + Z(tp,XcMB,0) 

Equation ( p3| ) shows that we can write 

1 + z(~ip,XcMB,0) = 1 + z (iP,Xcmb) + zi(iP,Xcmb) cos9 

where 
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1 + Zq(iP,Xcmb) 
Zi(iP,Xcmb) 



Rq 

Wo 



R r 



2a 
cA 



(1 - COS 1p COS Xcmb) 



2a R . 
cA Wq 



u(ljj) Rq 



A Wq 



Sin Xcmb 



(35) 



(using (|T^) in the last equality and assuming a < 0). The mean redshift of the CMB surface is zq; Z\ gives rise to 
an anisotropy in the CMB. We can therefore define the location of the CMB surface to be the t CMB (or Xcmb) that 
gives a mean redshift of 1000. That is, Xcmb is the solution of 



zq(4>,Xc 



1000 



(36) 



for any observer position ip. 

Having found Xcmb we can evaluate the anisotropy in the temperature of the CMB. Since T Q b s depends on the 
reciprocal of 1 + z the dipole moment in z will give rise to higher multipoles when expanded as a binomial series: 



Tobs(#) = 



1 + zo 



Zl 



l + z 



■ cos 9 



Zl 



1 + zq 



that is, 



8T{9) 



Zl 



1 + zo 



■ cos 9 



Zl 



cos 2 9 + 0(cos 3 



0(cos 3 9). 



(37) 



The dipole moment of the CMB temperature is then (using (p6|)) 



Si 



Zl 



l + za 



10" 3 zi. 



(38) 



Measurements of the CMB can now be used to constrain the model parameters. We at least require that the dipole 
moment should be no larger than the observed dipole anisotropy, \6i\ < Ti/Tq « 10 -3 (i.e., \z\\ < 1). If this is satisfied 
for any model then it is clear that the quadrupole and higher multipole moments will all be < 10~ 6 - certainly no 
larger than their observed values. In fact, such a constraint on zi is very weak, leaving vast tracts of parameter space 
entirely untouched. Moreover, there are very good reasons for believing that there is a significant contribution to 
the observed dipole moment from the peculiar velocity of the Local Group as a result of infall towards the Great 
Attractor [63, which can be measured with moderate accuracy using galaxy surveys, and seems to be consistent with 
the motion of the Local Group with respect to the CMB j6^|7(J, although see [7l] ]. Actually, it is not beyond the 
bounds of possibility that genuinely inhomogeneous background models such as those we consider h ere co uld mimic 
the grosser features of the local anisotropies (in particular the GA induced dipole effect - see Sees. |[IIA| and [fflj]) 
as well as the CMB dipole. The 'real' universe would then be a perturbation of this inhomogeneous background. 
Viewed in this way, local observations (galaxy surveys, etc.) would reveal the effects of the large-scale inhomogeneity 
and also contain information about smaller-scale perturbations (peculiar velocities and the density contrast). That is, 
the peculiar velocity field and the density contrast that are inferred assuming an FLRW background would actually 
be thought of as containing one part that reflects the difference between the inhomogenous and FLRW background 
models (and therefore contains the dipole effect referred to above, for example) and another part that is the 'true' 
perturbation. (The distinction between these two components is a fine one: ultimately it reflects the difference between 
linear and fully nonlinear perturbations of FLRW models.) We will consider this in a future paper j72|. However, 
we choose here to reject iconoclasm in favour of the more conservative viewpoint that most of the observed dipole is 
due to the peculiar motions induced by local inhomogeneities, but that there remains some leeway - up to 10% of the 
observed dipole - due to observational uncertainties, for there to be a purely cosmological contribution to the CMB 
dipole. Then the largest dipole moment that we can accept from our models is \Si\ < 10~ 4 , or 

12*1 < 0.1. (39) 

Given any model parameters and some observer position we adopt the following procedure. First we use ( p7|) to 
determine the epoch of observation for some Hq, as usual, then we solve for Xcmb using ([36]). Having found all the 
parameters we need to determine z\ we simply check ( 39[ ) to see whether the model, or at least that observer position, 
must be rejected. In practice we can simply solve (B9) as an equality to obtain a as a function of b at the boundary 
of the allowed region, and this is what is shown in Fig. |8| (for ip = tt/2), where it can be seen that a low value of Hq 
constrains our models most - in contrast to the age constraint. This is because Hq decreases monotonically with time, 
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FIG. 8. The exclusion diagram from CMB anisotropies for an observer at ij) = n/2. The curves for Ho = 50 and 
70 km s" 1 Mpc -1 are shown. As Ho increases, the 'fingers' move down to more negative a. The allowed region lies be- 
tween and above the fingers. 

so small Ho corresponds to a later time of observation and therefore a later time for the CMB surface, which means 
that W has evolved to become more inhomogeneous. We choose ip = ir/2 because, as is clear from ( |ll| ) and fl35|), the 
anisotropy is generally worst there, so if a model is rejected at ip = ir/2 it will be unacceptable everywhere. Again, 
therefore, for models not excluded in Fig. || detection of a CMB anisotropy of the magnitude that we observe would 
be typical, and the Copernican principle need not be abandoned for these models to be viable. 

The finger-like excluded regions in Fig. |^ appear because for different model parameters Xcmb takes on different 
values, and for some parameters this value is very close to (a multiple of) tt, so that the entire CMB seen by an 
observer in such a model is emitted from virtually a single point (the antipode). Since redshift depends only on the 
relative conformal factors at emitter and observer for our conformally flat models, the CMB must be almost exactly 
uniform however inhomogeneous the model (i.e., whatever the value of a). 

E. The Local Dipole Anisotropy 

Although we are not in a position to use real observations to constrain the dipoles that would be detected in 
observations of the 'local' universe (in galaxy surveys, for example, where z < 0.01, or with type la supernova data, 
for which z < 1), we can at least consider these effects qualitatively. It will turn out that for the models we consider 
observations of the dipole variation in the distance of objects at a given redshift (or the dipole in z at a given distance) 
directly constrain the acceleration of the fundamental observers. In fact, it is clear from ( |30| ) that this is the case for 
any model (at low redshift at least), so that measurements of these dipole moments at low z permit the acceleration u a 
to be measured. It is then only a question of the tightness of the constraints imposed by real observations, and the 
extent to which the acceleration signal can be separated out from other components in the peculiar velocity field. This 
we consider in [Q. Note that the acceleration dipole grows with distance at low z, as can be seen from (|30"|), so that 
it is distinct from the dipole resulting from GA infall, which amounts to an overall Lorentz boost by a constant factor 
and is therefore independent of distance (see Q). It should be noted that this is only true when objects more distant 
than the GA are observed; closer in the dipole structure due to GA infall is more complicated, and at very small 
distances there should be no dipole variation in iJ™ (to first order), because ii — in the standard interpretation 
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(although there will be a quadrupole component due to shear). 

In this section we adopt the null hypothesis that the universe is well described by a (perturbed) FLRW model 
(i.e., that observations are at least consistent with it = 0), and determine whether, and under what circumstances, 
local observations may provide a tighter constraint on the Dabrowski models than the CMB (Fig. ||). 

From (|2^) and (|l^) the redshift dipole for objects at any coordinate radius X from an observer at tp (in ( |l2| ) X is 
the observer position relative to the centre of symmetry, which we now denote ip) is seen to be 

\z 1 (iP, X )\ = ^P ^ |sinxi (40) 
c vAWb 



(cf. equation |35|). For small z this just corresponds to the H™ dipole in (p0|), since z = (Hq 1 / c)r prop (at low z the 
distance measures ta and tl are the same as proper distance rprop). It is clear that the redshift anisotropy for objects 
at some distance from the observer is closely related to the acceleration u of the fundamental observers and the dipole 
moment in the measured H™ in d30|). Defining SH = ii/c, the dipole in redshift for objects at a given x becomes 

\ziW,X)\ = — ^-\sm X \. (41) 

If we assume that the time of observation, T, is fairly close to t = (as is generally the case for the values of Hq 
we allow), then models with a < have T <C —bjla (the time at which the scale factor reaches its maximum value 
and H = 0) and Wo ~ 1, so that Rq/c w 6T and R,t /c ~ b. Equation ([0]) then becomes, with the help of (p7[), 

\zi(il>,X)\ « — 7=|sinx| « -7=1 sm X | — . (42) 
cyA no VA -Ho 

Note that dependence on distance from the observer only arises through the sin% factor, so for objects at any given 
distance from the observer 

\ Z1 W, X )\<JL*E. (43) 
v A t±Q 



In principle, therefore, local observations of redshift anisotropies determine the ifo-dipole, and therefore the acceler- 
ation u, through ( p0| ) and (|4^). Conversely, any measurement of the acceleration or Hubble dipole can be used to 
constrain the dipole anisotropy z\ at all redshifts through ( f42| ) or (|i^). 

Before con sideri ng the constraints imposed by local observations, we compare (|43|) with the CMB constraint (p9) im- 



posed in Sec. Ill E . If we assume that the bound in (p9|) is reached, i.e., that the CMB dipole due to the inhomogeneity 



of the Dabrowski model we are considering is \z\\ = 0.1, then it follows from ( J43| ) that 

(\f~K/b is a decreasing function of b, so that this constraint is stronger for smaller b). That is, if the universe was well 
described by a Dabrowski model and the cosmological CMB dipole was measured to be as large as the bound specified 
in d39|), then the local -f/o-dipole must be larger than O.l-v/A/6- (For example, a variation in Hq of less than 20% 
can only arise if vA/6 < 2, or b > l/s/5 « 0.45; for a variation of less than 10% we must have b > l/\/2 ~ 0.71.) 
Since the CMB anisotropy constraint (|3^) allows models with b > 0.5 (see Fig. |^), corresponding to a variation in 
Hq around the sky of at least 17% for models on the boundary of the allowed region in Fig. p where \zi\ = 0.1, if 
present observations show that the _ffo-dipole is less than 17% some of the models allowed by (39|) will be excluded 
by these local observations. Note, though, that (^) is considerably stronger than is really required for most model 
parameters, owing to the fact that the sinx factor in (|42 ) was neglected. This amounts to adopting as the CMB 
constraint the envelope of the fingers in Figs. ||, |^ and [L0|, which would obviously overlook large areas of parameter 
space that should really be allowed. 

Although it would seem that the easiest way to constrain the acceleration in the Dabrowski models would be to 
measure the local 77 -dipole, it is not really possible to measure SH/H accurately |72[ |. Moreover, there is known 
to be a dipole in observations of galaxies at somewhat larger redshifts (z ~ 0.01), which is usually interpreted as the 
effect of infall of the Local Group towards the Great Attractor |^,^9| . This infall is manifest as a systematic relative 
motion of the Local Group with respect to distant galaxies at a velocity v ~ 600 km s _1 , coinciding with the Local 
Group motion relative to the CMB frame. In order not to conflict with these observations we at least require that at 
redshifts zq w 0.01 the dipole moment due to the Dabrowski acceleration is no larger than the observed dipole: 
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czi < v — 600 km s 1 . 
Then, since at low redshift the linear Hubble law is valid, we have 

cz a = H r plop , czi = SHr prop , 

and so 

SH cs^ eOOkms-i 

H ~ cz ~ 0.01c ~ U 1 &j 

That is, we can allow at most a 20% variation in Hq around the sky. This is comparable to the 17% variation required 
by the CMB limit derived above, which means that (EjJ) does not provide a stronger constraint on the Dabrowski 



model parameters than (p9p. Interestingly, JT3 74 estimated the dipole in Hq to be about 20% (the Rubin- Ford effect). 
However, this measurement has since been discredited {75|, being the result of selection effects. 

If the _ff -dipole could actually be measured, or at least bounded, then the time independence of the acceleration 
scalar in (jL2|) makes it a simple matter to use this to constrain the model parameters a and b. From (^) and jl^ ) 
(with x = i } ) 

SH u 2\a\ . , 
tt = —rr = —f= — smip. 
H cH y^AHo 

If SH/Hq < 7, say, and we assume this holds for all observer positions tp, in accordance with the Copernican principle, 
then 

\a\ < i 7 ^oVA. 

It should be borne in mind that throughout the preceding discussion we have only considered the modulus of the 
dipole, not its direction. It is clear from (g5[) that the local dipoles may be in the same direction as the CMB dipole or 
in the opposite direction, depending on the sign of sinx CMB . What is more, the variation of the dipole with distance 
is controlled entirely by sinx, so the dipole will change sign whenever \ is a multiple of tt. The alignment of the 
CMB and local dipoles is usually taken to be a strong sign that both result from the peculiar motion of the Local 
Group (due largely to GA infall). However, it is clear that for the models we consider these dipoles will also always 
be aligned (or anti-aligned). 



F. The Combined Exclusion Diagrams 

When we combine all of the constraints derived in this section (Figs. || and |l^) we can see that for Hq = 
50 km s _1 Mpc -1 the strongest constraint comes from the CMB, with age placing somewhat weaker limits on the 



allowa ble de gree of inhomogeneity (which is measured largely by the size of a - see Sec. IV). The 'size' restriction 



of Sec. Ill C eliminates quite a large region of parameter space for small b, but this is not really a constraint on the 
inhomogeneity, which is our principal concern. 

Perhaps rather surprisingly, given the results of [ [32"|]35| ], the strongest constraint for larger Hq really comes from 
the age. As can be seen from the exclusion plot for Hq = 70 km s _1 Mpc -1 , Fig. |l^, the CMB constraint pokes 
out in places to eliminate certain regions, and the size constraint cuts off low values of b, but age does most of the 
dirty work. It can also be seen that for Hq = 80 km s _1 Mpc -1 age imposes a very strong constraint on the models 
(dashed line in Fig. |l0| ): models with high Ho must be very nearly homogeneous, or they are too young. However, if 
we relaxed the age constraint to 10 Gyr the CMB anisotropy would be the dominant limitation for most values of Hq 
in the currently fashionable range (50 < Ho ^ 80 km s" 1 Mpc -1 ). 

We should not forget, at this point, to reintroduce the restriction (pp| ) from the dominant energy condition, which 
rules out high b. This is not shown on the diagrams, in order to avoid clutter. Most of the models eliminated by this 
constraint have already been ruled out by the age or CMB constraints, and models that are rejected solely by the 
dominant energy condition are not hugely inhomogeneous (see the next section). 
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b 




FIG. 9. The complete exclusion diagram for all the observational constraints studied (age, size and the CMB anisotropy), 
for Ho — 50 km s _1 Mpc" 1 . We have taken the 12 Gyr age constraint, to be conservative. The dominant energy condition 
should be added to these constraints: it eliminates models with b > 0.82 (equation (ppf)). 

b 



0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 




FIG. 10. As in Fig. but for Ho = 70 km s 1 Mpc 1 . The age constraint for Ho = 80 km s 1 Mpc 1 is also shown as a 
dashed line (close to the 6-axis): high Ho means a low age. 



20 



IV. THE SIZE OF THE INHOMOGENEITY 



While we have considered many different aspects of the Dabrowski models, what we have not done is to assess 
the extent to which the models that are not excluded are inhomogeneous. It is obvious from the exclusion plots, 
Figs. |§ and [To], that the homogeneous Dabrowski models (those with a = 0) are the 'most acceptable', in that all the 
constraints favour small a. This should not be surprising, as far as anisotropy constraints are concerned, at least. 
What is not clear is whether the allowed region only contains models that are very nearly homogeneous. We will show 
that it does not. 

The most natural way to assess the degree of inhomogeneity of the models is to examine the variation of the 
pressure over surfaces of constant cosmic time. It can be seen from ( p"7| ) that the extremes of pressure occur at x = 
and x = 7r, so we define the inhomogeneity factor II to be the relative pressure difference between the two poles: 



n = 



p(n) -p(0) 



p(0) 



W. (46) 



If IT > 1 then it is reasonable to say that the models are truly inhomogeneous, whereas if II <C 1 they are obviously 
nearly FLRW. Note, though, that II depends on the cosmic time surface under consideration: for small t, H ss 0, 
and II reaches its maximum at t = —bjla (see (H)), at which time 

26 2 1 - A 

n =I-^ = 2 ^ 

so that the models are significantly inhomogeneous (II > 1) when 

6 > -i= « 0.58. (47) 

Most of the allowed mode ls in Figs. [)] and [l(] satisfy ( f47j ) - models with smaller b have already been eliminated by 
the size constraint in Sec. Ill C| . 



This is not really a fair reflection of the inhomogeneity of the models at the times of observation that are relevant 
here, though, because the Hq constraint d27j ) generally ensures that the epoch of observation is quite early on in the 
evolution of the universe when the scale factor is somewhat smaller than its maximum size. To evaluate the impact 
of this, consider two specific examples. From the allowed region of Fig. O choose the model at a = —7, b = 0.75. 
For H = 50 km s" 1 Mpc" 1 the solution of @ is T sa 0.016 Mpc s km~^ which gives II = 1.33. For the model at 
a = -8, b = 0.64 we get T = (3 - \/5)/50 w 0.015 Mpc s km" 1 and II = 0.86, which is close enough to 1. These 
models are certainly not 'close to FLRW. 

Nevertheless, it could be said that the models are not massively inhomogeneous, and that their degree of in- 
homogeneity only reflects the looseness of the constraints applied. This is not so: even for models that are only 
inhomogeneous at the 10% level (II = 0.1), the CMB anisotropy, Si < 10~ 4 , is at least three orders of magnitude 
smaller than the inhomogeneity that generates it. This certainly conflicts with the spirit of the almost EGS theorem 
of |^| which says that small CMB anisotropies indicate correspondingly small perturbations from homogeneity. Of 
course, we know from JIJ that the Dabrowski models all admit a perfectly isotropic radiation field; the fact that the 
CMB is not precisely isotropic here is a result of our more sophisticated definition of the CMB surface. 



V. A NOTE ON MODELS WITH a > 



In this paper we have only considered the effects of the models with a < in order to satisfy the strong energy 
condition. The type la supernovae data strongly implies that the expansion rate of the universe is increasing (go < 0) 
which, in the A = FLRW case at least, can only occur when the strong energy condition fails.as may be seen from 
the Raychaudhuri equation. This is not necessarily true when acceleration is present however [|6||8j] , but we may derive 
the following condition for the Dabrowski models to have go < at some time t: 



A + ((2at + bf + b 2 ) sin 2 ~ 



> 0. 



The term in brackets is always positive (unless A < 0, which we do not consider), so that q < if and only if a > 0. 
In fact, this is intuitively clear from the quadratic form of R(t): the accelerated expansion implied by the SNIa data 
can only be produced when R(t) is an 'upright' quadratic (a > 0). 
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By analogy with FLRW models, we may define a Hubble normalised density parameter 



which we may use to get a rough idea of how these models behave in comparison with FLRW models. The present 
day value of the density parameter is roughly in the range 0.3 < Oq < 1. For the Dabrowski models, we note that in 
the limit t — > we have SI — > 1/b 2 , and H t t — > — 4a/6 3 . We know that b is positive in order that the scale factor is 
positive, which implies that fi^ is negative at the big bang only if a > 0. However, we also see that as t — > 0, then 

approaches some value larger than 1 if b < 1 (which is required if we reject models with finite-density singularities). 
This means that if we desire f2 ^ 1 today, then we need a > 0, in contradiction with the strong energy condition. 

When a > 0, the quantitative restrictions from the age, size and CMB anisotropics are not strongly affected; 
indeed, in some ways the outlook is better. As is shown in the proper age of an observer in an a > model is 

considerably larger than their compatriot in a negative a model, for the same value of Hq; the 'size' of the universe 
(i.e., the redshift of the antipode) is also larger for a > 0. The CMB anisotropy goes as \z\\ ~ \a\, which suggests that 
it is unaffected by the sign of a.[] 

However, a > implies that the universe will open up at t = (1 — b)/2a, after which time the universe is infinite in 
spatial extent (note that models with a > expand forever). Spatial infinity occurs at some finite ip = ip malcl where 
W = in (|l^). Applying the Copernican principle in this case becomes difficult - at least in terms of producing 
exclusion diagrams. We will not give further attention to this here. 



VI. CONCLUSIONS 



We have studied the physical, geometrical and observational characteristics of a two-parameter family of inhomo- 
geneous, perfect fluid cosmological models which form a subclass of the models admitting isotropic radiation fields 
found in JjJ and coincide with the model I spacetimes of |54j. We have shown that these models do not suffer from 
particle horizons. 

The inhomogeneity of these Dabrowski models makes the investigation of their observational characteristics from 
any position essential, and the simple conformal geometry of the models was used in Sees. [n HI to derive exact, 
analytic expressions for the observational relations and other properties for any observer position; the constraints 
imposed by observations were examined. For any model parameters and any observer position we fixed the coordinate 
time T of the present epoch (i.e., the time of observation) by demanding that the value of Hq at that T, given by J27|), 
was consistent with current estimates of the Hubble constant ( p8|) . We then proceeded to test the consistency of the 
models with a variety of observational constraints, the most important being age and the anisotropy of the CMB. 

Obviously, in any viable cosmological model the age of the universe must be greater than the estimated age of of 
any of its components by the time the expansion has slowed to the rate we measure today. This is position dependent 
in our models: on a given surface of constant coordinate time (determined by the expansion rate Hq) the observers 
at the centre will generally be older than their counterparts elsewhere (see Fig. [l]). As with all the tests we apply, we 
show the constraints for the 'worst case' observer position, i.e., that which excludes the largest region of parameter 
space. 

Despite the fact that the Dabrowski models admit an isotropic radiation field H , this will not in general correspond 
exactly to the actual cosmic background radiation field after decoupling. This is because, although the models are 
homogeneous at early times (conformal factor W = 1), the inhomogeneity (W — W(x,t)) that develops up to the 
time of decoupling will leave its imprint in the CMB. In other words, the EGS-type theorems prove the existence of an 
isotropic radiation field, but whether this is realised in a particular spacctimc depends on the physics of decoupling. 
We used a pragmatic definition of the time of decoupling to estimate the effect of the inhomogeneity on the CMB. 
The anisotropies in the observed CMB temperature go as 1/(1 + z) ~ 1/W: that is, all anisotropies arise entirely from 
the inhomogeneities in the conformal factor W at the time of decoupling; the conformal flatness of the spacetimes 
allows the light from the CMB surface to travel unmolested by any subsequent deviation from homogeneity, and we 
see the CMB as a reflection of the universe's inhomogeneity at that time. Since t is small at decoupling W ~ 1 and 
the inhomogeneities are very small. The dipole moment of the CMB is the largest moment, which we restrict to be 



1 This is not quite true, however, because it depends on the time at which the surface of last scattering occurs (i.e., on Xcmb)- 
This will be later for models with a > (since the epoch of observation is later, for the same Ho), and the universe will be 
more inhomogeneous. However, this effect turns out to be rather weak and does not affect the conclusions we reach. 
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smaller than 10% of the observed dipole, on the basis that present data cannot show that Local Group motion and 
the CMB dipole agree to within 10% @. 

In addition to these observational constraints, we required that the first zero of the distance-redshift relations does 
not occur too close to the observer, on the grounds that this would probably have been observed. This ruled out 
models with small b, which turned out to be fairly homogeneous anyway. 

The matter content of these models is unusual, in that there is no equation of state. However, we know from |IJ (see 
also ||,[58]|) that they do admit a thermodynamic interpretation (i.e., definitions of number density, temperature and 
entropy throughout spacetime consistent with the first law of thermodynamics). We have shown that many, though 
not all, of the Dabrowski models satisfy the weak, strong and dominant energy conditions (Sec. II D), so that their 
matter content cannot be ruled out as obviously unphysical. We restricted our detailed analysis to the models that 
satisfy all three energy conditions (the strong energy condition then requires a < 0). Of course, we are not forced to 
accept the energy conditions. Although they seem physically very reasonable conditions to impose on any form of 
matter there are many examples in cosmology of matter that does not satisfy them (quantum fields, for example, can 
exhibit negative energy density). Given that quintessence models ]44|-[47|] , which produce an accelerating expansion 
in FLRW models (consistent with the supernova data), must break the strong energy condition, rejection of the 
a > models may be premature, and we have considered their properties briefly in Sec. ^ (see also ||). 

Our studies have shown that there is a significant subset of the Dabrowski models that are markedly inhomogeneous 
but cannot be excluded on the basis of the tests considered here. It is possible, for every observer in each of the 
models in the allowed regions of Figs, [s] and [l(], to choose the epoch of observation so that the observed value of H 
is reproduced, the age is greater than the measured age of the universe and there are no obviously unacceptable 
features at low redshift (z < 5) in the observational relations. Most importantly, though, the dipole in the cosmic 
background radiation would be considerably smaller than the observed CMB dipole: despite the inhomogeneity of the 
models the anisotropy they produce is very small. The fact that this is true for every observer means that it is not 
possible to reject these models by appealing to the Copernican principle. As a result, the standard assumption that 
the observed high degree of isotropy about us combined with the Copernican principle necessarily forces the universe 
to be homogeneous (i.e., the cosmological principle) is seriously undermined. 
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APPENDIX A: TRANSFORMATION TO A NON-CENTRAL POSITION 

We want to transform the metric (jj^) from the (x, 9, (j>) coordinate system, whose origin is at the centre, to coordinates 
centred instead on some observer at \ = while preserving the form of the FLRW part of the metric. It is therefore 
necessary to identify the transformations of the (homogeneous) FLRW spatial sections that leave the FLRW metric 
invariant, i.e., the isometries of the spatial sections. This is simple. Since the spatial sections of an FLRW model with 
positive curvature constant (A > 0) are 3-spheres, the isometries we require are 4-dimensional rotations (i.e., elements 
of 5*0(4), the isometry group of the 3-sphere). 

A sphere of radius R in 4-dimensional space with cartesian coordinates (x, y, z, u) is defined by 

x 2 + y 2 + z 2 + u 2 = R 2 . 

We have three coordinates on this sphere: x an d the two spherical polar angles 9 and <fr. These are related to the 
cartesian coordinates by 

x = R sin x sin 9 cos tj> (Al) 

y = R sin x sin 9 sin <fi ( A2) 

z = i? sin x cos # (A3) 

u = Rcos x (A4) 

The origin, x = 0, is then at x — y = z = 0, u = R. We are only interested in 4-rotations that move the origin, and, 
as the initial metric is spherically symmetric (really spherically symmetric, not just conformally: even the conformal 
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factor is spherically symmetric about the centre), we need only consider moving the observer in one direction, which 
we choose to be the z direction (i.e., to a position with non-zero z, but x = y = 0). Clearly, then, we are looking for 
a rotation in the u — z plane. Since we have the conformal factor as a function of x we want to find x as a function 
of the new coordinates. Starting with coordinates x', 0' and 0', centred on some position x = along with their 
prim e d ca rtesian counterparts x 1 ', y' , z' and u' (which are related in the same way as the unprimed coordinates in 



(Al)-(A4)), a rotation back to the original coordinates is given, in cartesian coordinates, by x — x' , y — y' and 

z = cos ijjz' + sinipu' , 

u = — sin-0z' + cos tpu'. (A5) 
(Note that at the origin of the primed coordinates, where z' = and v! — -R, we have u = Rcosip, showing that x = V- 1 



there, as required.) Equation (A5), along with the primed versions of (A3) and (A4), then immediately gives 



cos x = cos ip cos x' — sin ip sin x' cos 9' , ( A6) 

and this is all we will need, since the only spatial coordinate that enters into the original metric (0) is x, and that 
enters only as cos x (2 sin 2 ^ = 1 — cos x) ■ 

APPENDIX B: REDSHIFT IN CONFORMALLY RELATED SPACETIMES 

We present two derivations of the redshift formula for spacetimes sharing some of the simple properties of the 
Dabrowski models. The first can be used in spacetimes that are conformal to simpler metrics for which the geodesies 
and redshifts can be found. The second is valid for spacetimes that are conformal to a spherically symmetric spacetime. 

If we have two conformally related metrics, g a b = Q 2 gab (f2 > 0), their associated metric connections are related by 
(see appendix D of fiofl ) 

V b V a = V b V a + (V b In SI) V a + (V C V C In SI) 5 a b - (V d In SI) g ad g bc V c , (Bl) 

for any vector field V a . Null geodesies with respect to g (or, more correctly, with respect to V) satisfy g a bk a k b = 
and 

k b V b k a = 0, (B2) 
where k a is the tangent vector to the geodesic. Applying (p3l|) gives 



= k b S7 b k a = k b \7 b k a + 2k b \7 b lnSlk a -g ad S7 d \nSl (g bc k b k c ) = -^(Sl 2 k b )V b (Sl 2 k a ), 
and we see immediately that 

k a = Sl 2 k a (B3) 

is the tangent vector to a null geodesic with respect to V. So, every null geodesic of V corresponds to a null geodesic 
of V (and vice versa, since we can repeat the above steps interchanging g and g and putting SI t— ► I /SI). If the 
geodesies with respect to V are known explicitly then we can find them easily for V. 

To find the redshift, though, we also need the velocities of emitter and observer. A four-velocity satisfies g ab u a u b = 
— c 2 , and if we define 

u a = Slu a (B4) 
then u is a four-velocity with respect to g: g ab u a u b — —c 2 . Redshift is calculated from the ratio of the emitted 



frequency v E = u a k a \ E to the observed frequency v Q — u a k a \ Q . Using (B3) and (B4) we have 

u a k b 1 



which means that 



u a k a = g a bu a k b = Sl 2 g ab ^-^-r = ^-u a k a , 

Si Si Si 



u b k b \ Q Sl E u b k b \ o Sl E 
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where z is the redshift associated with g a b for the fundamental velocity u. If the paths of null rays in the spacetime g a b 
and the redshift formula for the velocity u are known then (B6) gives the redshift in the true spacetime. For the 
Dabrowski models g a b is an FLRW metric, the conformal factor is fl = 1/W (see (P)) and u is the usual FLRW 
comoving velocity field. The well-known expression for redshift in FLRW spacetimes, 1 + z = R /Re, then gives 



l + z = 



Ro_We 
Wo R E 



(B6) 



When the true spacetime is conformal to a spherically symmetric spacetime the radial null geodesies connecting any 
point with an observer at the centre are obviously purely radial (since their paths are not affected by the conformal 
factor). They are therefore given (in terms of coordinates r and t with respect to which the spherical symmetry is 
manifest) by some function t (r E , t E ) relating the time, t , that the light ray is received by the observer, to the time 
of emission, t E , for an object at radius r E . This is just the lookback-time relation. Redshift, as the ratio of proper 
time intervals dr at the observer to proper time intervals dr E at the emitter, is then given by 



j _ dT dr dt dr ( dt Q dr E dt Q dt E 
dr E dt dr E dt Q V dr E dr E dt E dr E 



1 



dt Q r dt 

o t 

4, \o^ Ue dT v 



(B7) 



(When the coordinates r and t are comoving - u r = - the r-derivative term disappears.) This will provide an analytic 
expression for the redshift whenever the lookback-time equation can be integrated. For the Dabrowski models: 



0. 



1/2 



W 



(B8) 



and the lookback time can be derived directly from the metric: on the past null cone of the observer ds = = dO = 
leading to an expression for dx/dt, which, when integrated, gives 



X 



nr f T \ dt' 



Differentiating this with respect to t at fixed \ then gives 

dto^dT = Rc 
dt E dt R E 



which, together with (B7) and (B8), results in the expression (B6) for the redshift 
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